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We present a comprehensive study of polar and magnetic excitations in BiFeOs ceramics and 
a thin film epitaxially grown on an orthorhombic (110) TbScOs substrate. Infrared reflectivity 
spectroscopy was performed at temperatures from 5 to 900 K for the ceramics and below room 
temperature for the thin film. All 13 polar phonons allowed by the factor-group analysis were 
observed in the ceramic samples. The thin-film spectra revealed 12 phonon modes only and an 
additional weak excitation, probably of spin origin. On heating towards the ferroelectric phase 
transition near HOOK, some phonons soften, leading to an increase in the static permittivity. In 
the ceramics, terahertz transmission spectra show five low-energy magnetic excitations including 
two which were not previously known to be infrared active; at 5K, their frequencies are 53 and 
56cm“^. Heating induces softening of all magnetic modes. At a temperature of 5K, applying an 
external magnetic field of up to 7 T irreversibly alters the intensities of some of these modes. The 
frequencies of the observed spin excitations provide support for the recently developed complex 
model of magnetic interactions in BiFeOa (R.S. Fishman, Phys. Rev. B 87, 224419 (2013)). The 
simultaneous infrared and Raman activity of the spin excitations is consistent with their assignment 
to electromagnons. 

PACS numbers: 78.30.-j; 63.20.-e; 75.30.Ds 


I. INTRODUCTION 

Among novel materials, an intense effort is concen¬ 
trated on the study of multiferroicsjii^ A wide range of 
applications, such as information storage, sensing, actu¬ 
ation, and spintronics, await pioneering materials and 
strategies that would produce robust magnetoelectric 
coupling at room temperature (RT)i^i^ The ability to ma¬ 
nipulate magnetization in a magnetoelectric multiferroic 
by electric fields can be extremely promising for such ap¬ 
plications, due to the simplicity and cost-efficiency of ap¬ 
plying an electric field. As one of the few single-phase RT 
magnetoelectric multiferroics, bismuth ferrite BiFeOs, is 
at the center of attention, as it presents a ferroelectric 
phase transition at approximately HOOK and an antifer¬ 
romagnetic one at 643 

The knowledge of lattice and spin excitations in 
BiFeOs is essential for understanding the underlying 
mechanisms that induce its multiferroic behavior. A se¬ 
ries of Raman and Infrared (IR) spectroscopy studies 
have presented controversial results concerning the as¬ 
signment of the magnon and phonon modes, as well as of 
the highly acclaimed electromagnons (i.e., electrically ac¬ 
tive magnons). Probing IR-active low-energy excitations 
is hindered by a lack of sufficiently large single crystals. 
Raman-active phonons of the rhombohedral i?3c BiFeOs 
structure have been reported for single crystals^^— poly¬ 
crystalline ceramics^^— and thin fllms )^^d^ however. 


there is a significant discrepancy between the claimed 
phonon frequencies and symmetry representations. A 
possible explanation for such inconsistences between the 
various experimental results is the presence of oblique 
phonon modes, which show a continuous variation of fre¬ 
quency along the phonon propagation vector with respect 
to the crystallographic axis of the probed specimeniii Up 
to now, the phonon IR spectroscopy studies have been fo¬ 
cused on ceramicsi^^— and single crystalswhereas, 
to our knowledge, no report of a thin-film IR investiga¬ 
tion exists. 

When it comes to magnon and electromagnon stud¬ 
ies, Raman spectroscopy holds the record for the num¬ 
ber of excitations observed, yielding up to 16 mag¬ 
netic modes corresponding to the cycloidal spin struc¬ 
ture (8 cyclon and 8 extra-cyclon modes) for sin¬ 
gle crystals but only one for polycrystalline thin 
filmsFurthermore, the influence of strain on the 
number of spin excitations and their frequencies were 
reported for epitaxial thin films grown on a series of 
substratesi^ Time-domain THz spectroscopy^ inelas¬ 
tic neutron scattering measurements21“— and absorbance 
spectroscopy in the THz ranged have revealed modes 
similar to the Raman-active onespredicted by theoret¬ 
ical calculationsThe best agreement between the 
experimental and theoretical spin-excitation frequencies 
(including their splitting in external magnetic field) was 
obtained for a microscopic model that takes into ac- 
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count the nearest and next-nearest neighbor exchange 
interactions, two Dzyaloshinskii-Moriya interactions and 
an easy-axis anisotropy^^*^ The same model success¬ 
fully described the low-energy inelastic neutron scatter¬ 
ing spectral^ In contrast, Komandin et ali^ showed IR 
transmission spectra with an excitation at approximately 
47 cm“^ which had not been previously reported by ex¬ 
perimental or theoretical studies. The intensive discus¬ 
sion in the literature concerning the nature of the spin ex¬ 
citations raises the question: which excitations are pure 
magnons (i.e., contribute only to the magnetic perme¬ 
ability /i) and which are electromagnons (i.e., influence 
at least partially the permittivity e*)? It is worth noting 
that according to a recent symmetry analysis,— BiFeOa 
allows directional dichroism and therefore spin waves can 
be simultaneously excited by the electric and magnetic 
components of electromagnetic radiation. For more de¬ 
tails on BiFeOa spin dynamics, see the review of Park et 
al^ 

In the current work, we report spin and lattice exci¬ 
tations in BiFeOa ceramics, as measured by the com¬ 
bination of IR reflectivity and time-domain THz trans¬ 
mission spectroscopy, in a temperature range from 10 to 
900 K. All 13 IR-active phonon modes are observed, ex¬ 
hibiting softening on heating. Five low-frequency spin 
modes are detected from 5 K up to RT, the highest two 
appearing at 53 and 56cm~^. This corresponds to the 
frequency range where such excitations were theoretically 
predictedbut not experimentally confirmed up to 
now. At 5 K, the low-energy spin dynamics in the THz 
range were also studied in a varying magnetic field of up 
to 7 T. Softening of the (electro)magnon frequencies upon 
increasing the magnetic field was observed. Additionally, 
a BiFeOa epitaxial thin film grown on an orthorhombic 
(110) TbScOa single crystal substrate was studied for the 
first time via IR reflectance spectroscopy. 


II. EXPERIMENTAL DETAILS 

BiFeOa ceramics were prepared by the solid-state 
route. A stoichiometric mixture of Fe 2 03 and Bi 2 03 
powder oxides with a purity of 99.99% was ground and 
uniaxially cold-pressed under 20-30 MPa pressure into 
8 mm diameter pellets. The pellets were then covered 
by sacrificial BFO powder to avoid bismuth oxide loss 
from the pellet, and sintered in a tube furnace at 825 °C 
for8h in air. To avoid any secondary phase formation, 
the samples were quenched to RT. Polished disks with a 
diameter of 4 mm and thicknesses of ca. 600 and 338 /xm 
were used for the IR reflectivity and THz transmission 
measurements, respectively. 

An epitaxial BiFe 03 thin film with a thickness of 
300 nm was grown by reactive molecular-beam epitaxy 
on a (110) single crystal substrate. The growth pa¬ 
rameters were the same as for the samples reported in 
Ref. [H. Near-normal incidence IR reflectivity spectra 
of the BiFe 03 ceramics and film were obtained using 


a Fourier-transform IR spectrometer Bruker IFS 113v 
in the frequency range 20-3000 cm~^ (0.6-90THz) at 
RT; for the low and high temperature measurements the 
spectral range was reduced to 650 cm“^. Pyroelectric 
deuterated triglycine sulfate detectors were used for the 
room- and high-temperature measurements up to 900 K, 
whereas a He-cooled (operating temperature 1.6 K) Si 
bolometer was used for the low-temperature measure¬ 
ments down to 10 K. A commercial high-temperature cell 
(SPECAC P/N 5850) was used for the high-temperature 
experiments. The thermal radiation from the hot sample 
entering the interferometer was taken into account in our 
spectra evaluation. 

THz measurements from 3cm“^ to 60cm“^ (0.09- 
1.8 THz) were performed in the transmission mode 
with the use of a custom-made time-domain terahertz 
spectrometer. In this spectrometer, a femtosecond 
Tiisapphire laser oscillator (Coherent, Mira) produces 
a train of femtosecond pulses which generates linearly 
polarized broadband THz pulses in a photoconduct¬ 
ing switch TeraSED (Giga-Optics). A gated detection 
scheme based on electrooptic sampling with a 1 mm 
thick [110] ZnTe crystal as a sensor allows us to mea¬ 
sure the time profile of the electric field of the trans¬ 
mitted THz pulse. The same high-temperature cell as 
for the IR reflectivity was used for the THz range high- 
temperature measurements. Oxford Instruments Opti- 
stat optical cryostats with mylar and polyethylene win¬ 
dows were used for the low-temperature THz and IR mea¬ 
surements, respectively. THz experiments in an external 
magnetic field i7ext < 7 T were performed upon decreas¬ 
ing H with an Oxford Instruments Spectromag cryostat 
in the Voigt configuration, where the electric component 
of the THz radiation Ethz was set perpendicular to i7ext. 

The IR reflectivity and THz complex (relative) per¬ 
mittivity spectra £*(w) were carefully fit assuming the 
factorized form of the dielectric function based on a gen¬ 
eralized damped-harmonic-oscillator models^ 


e*(w) = £oo n 


N 2 
^LOj 


‘^TOj 


-I- iwynoj 
-I- iwyTOj ’ 


( 1 ) 


where wtoj £md wloj the frequencies of the j-th 
transverse optical (TO) and longitudinal optical (LO) 
phonons, qxoj and ynoj are the corresponding damp¬ 
ing constants, and Eoo denotes the high-frequency (elec¬ 
tronic) contribution to the permittivity, determined from 
the RT frequency-independent reflectivity tail above the 
phonon frequencies. The reflectivity i?(w) is related to 
the complex dielectric function e*{u}) by: 


R{uj) 


V£*i^) - 1 


( 2 ) 


To evaluate the IR reflectance spectra of the 
BiEe 03 /TbSc 03 thin film, a model corresponding to a 
two-layer optical system was used4^ The IR reflectivity 
spectra of the bare TbScOs substrate were fit first at each 
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FIG. 1: (Color online) (a) IR reflectivity spectra of BiFeOs ce¬ 
ramics at selected temperatures; (b) Real, and (c) imaginary 
part of permittivity, as obtained from the fits. 


temperature, using the model described by Eq. [TJ The 
values of oscillator parameters (wTOjj 7TOjj 7LOji 

and £oo) obtained for the substrate were used for fitting 
the BiFeOa/TbScOa thin-film spectra. The model used 
to describe the complex permittivity of the thin film is 
a sum of N independent three-parameter damped har¬ 
monic oscillators, which can be expressed asi^ 


e*(a;) — Soo + 


N 


Aej-w^oj 


H + tW7TOj ' 


( 3 ) 


where Asj is the dielectric strength of the j-th mode. 
Eq. [3] is simpler than Eq. [U and its use is well justified 
because the damping of the LO phonons of the film does 
not appreciably influence the reflectance spectra. 


III. RESULTS AND DISCUSSION 

A. Phonons in the BiFeOs ceramics 

The temperature dependence of the experimental IR 
reflectivity spectra of the BiFeOa ceramics is shown 


in Fig. [TJa). The reflectivity band intensities are 
higher than those in the previously published spectra 
of ceramicsi^^— and single crystal,— but comparable to 
the single-crystal spectra published by Lobo et al— This 
confirms the very high quality and density of our ceramic 
samples, which is essential for an accurate determination 
of phonon and magnon parameters. 

The IR reflectivity and THz transmission spectra were 
fit simultaneously using Eq. [T] the resulting phonon pa¬ 
rameters are listed in Table 1. The real and imaginary 
parts of the complex dielectric spectra calculated from 
the fits are shown in Figs. [TKb),(c). As predicted by the 
factor group analysis for the rhombohedral i?3c structure 
of BiFeOaall 13 IR-active phonons (4Ai + 9E sym¬ 
metries) are clearly seen from 10 K up to RT. In Table 1, 
the symmetries of all modes are assigned based on Raman 
spectra^ first-principles calculations?^ and IR spectra 
of single crystals— The damping of the modes strongly 
increases on heating; therefore, above 300 K, the IR re¬ 
flection bands broaden and mutually overlap. Note the 
remarkable softening of the modes below 150 cm“^ upon 
heating (see Fig. [T]) . Consequently, the static permittiv¬ 
ity e(0) determined by the sum of phonon contributions 
as £(0) = £oo + X) increases on heating in agreement 
with the Lyddane-Sachs-Teller relation^ (see Fig. A 
qualitatively similar temperature dependence was pub¬ 
lished in previous studies of BiFeOa ceramicS )^^d^ but 
the absolute value of permittivity was lower, probably 
due to a lower density/quality of the earlier studied ce¬ 
ramics. In contrast, the permittivity calculated from the 
single-crystal spectra— is slightly higher, due to the ab¬ 
sence of the Ai modes whose dielectric strength Asj is 
lower than the F-symmetry modes. 

We would like to stress that the temperature depen¬ 
dence of permittivity calculated from the phonon con¬ 
tributions has been previously published mostly below 
rT ^i5,i 7 pjjjy reported £(0,T) up to 900 K. 

Experimental low-frequency (i.e., below 1 MHz) dielec¬ 
tric data are also not available above RT due to the pres¬ 
ence of significant leakage currents in BiFeOa at high 
temperatures. There is only one publicationi^ present¬ 
ing IR spectra of BiFeOa ceramics up to 1280K, i.e., 
far above Tq, however, in this case, emissivity was mea¬ 
sured instead of reflectivity above 600 K and £(0,T) was 
not reported. In that work, Massa et alt^ listed phonon 
parameters at selected temperatures, from which we re¬ 
calculated £(0,r); these values are presented for sake of 
comparison in Fig. [21 One can see that, at T = 523 and 
850 K, Massa et alt^ obtained the lowest published val¬ 
ues of £(0). We also note the unusual increase in £00 and 
£(0,T) on cooling below RT; the latter appears to be in 
contradiction with the soft-mode frequency increase on 
cooling, reported in the same work. Finally, two addi¬ 
tional phonon modes at T = 4K, not predicted by the 
factor-group analysis, are reported in Ref. [H. All these 
facts raise questions about whether these results describe 
the intrinsic properties of BiFeOa. Nevertheless, it is 
worth noting the peak in Massa’s data at the ferroelec- 
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TABLE I: A comparison of the parameters of the IR-active modes in the BiFeOa ceramics (at 10 and 900 K) and the thin film 
(at 10 K), obtained from IR spectra htting. The symmetry assignment of each mode is also given. 


Ceramics Thin film 


Symmetry 


10 K 



900 K 


10 K (Eext II [110]) 

10 K (Eext II [001]) 

^ CJTO 

7TO 

^LO 

7LO 

CJTO 

7TO 

CJLO 

7LO 

CJTO 

7TO 

Ae 

CJTO 

7TO 

Ae 

E 

75.8 

4.6 

81.5 

2.8 

64.3 

40.9 

70.5 

19.3 

74.7 

1.4 

4.3 

77.8 

1.8 

4.5 

E 

133.7 

22.8 

137.4 

34.6 

101.0 

65.1 

117.9 

51.1 

132.2 

2.3 

8.6 

136.7 

2.9 

12.2 

Ai 

137.8 

13.6 

175.8 

4.0 

124.3 

57.6 

165.8 

32.9 

145.2 

2.4 

7.4 




Ai 

231.2 

14.0 

233.9 

5.6 

226.2 

59.6 

233.5 

31.4 

225.9 

2.9 

0.6 

227.4 

2.2 

0.7 

E 

236.7 

7.6 

243.6 

13.3 

234.0 

31.7 

240.8 

69.3 

241.0 

4.3 

0.9 

242.4 

4 

0.2 

E 

264.0 

13.0 

283.0 

38.0 

258.2 

60.4 

278.0 

48.5 

266.2 

2.3 

5.4 

266.8 

4.4 

0.4 

E 

288.5 

26.2 

345.8 

9.2 

278.4 

61.6 

340.7 

40.4 

293.9 

4.6 

3.6 

279.5 

5.5 

10.6 

E 

349.0 

10.0 

368.1 

5.0 

342.5 

42.5 

363.6 

41.0 







E 

372.0 

7.3 

430.7 

10.6 

364.5 

50.1 

420.3 

51.4 




371.7 

9.2 

0.9 

E 

439.6 

11.0 

468.0 

13.3 

424.0 

54.3 

464.3 

37.1 




445.7 

0.3 

1.1 

Ai 

476.2 

41.2 

503.7 

49.1 

469.8 

41.0 

502.8 

106.2 

472.2 

7.1 

0.01 




E 

520.2 

21.8 

523.8 

76.0 

514.7 

45.7 

518.3 

62.5 

524.0 

9.8 

0.3 

516.8 

12.7 

0.04 

Ai 

549.5 

36.9 

606.2 

27.1 

542.1 

78.8 

595.3 

102.3 

550.6 

7.2 

0.07 

551.8 

2.3 

1 

Spin wave 











603.6 

9.9 

0.01 



Temperature (K) 

FIG. 2: (Color online) Temperature dependence of the static 
permittivity calculated from the sum of phonon contributions. 
The data for the ceramics and the film (Eext || [110]) from the 
current work are compared with previously published data 
obtained, in order, on ceramics,— and a single crystal.— 
The line connecting the data from Ref. [l^ is a guide for the 
eyes. See the text for additional details. 


trie phase transition near 1120 K; this is the only experi¬ 
ment showing a dielectric anomaly at Tc in BiFeOs. The 
peak is much lower than in other canonical ferroelectrics, 
apparently due to the first-order character of the phase 
transition from i?3c symmetry to the high-temperature 
Pnma structurei^i^ 


B. Excitations in the BiFeOa thin film 

The polarized IR reflectivity spectra of the bare 
(110) TbScOa substrate and reflectance spectra of the 
BiFeOa/TbScOa thin film are shown in Fig. All 
BiFeOa mode frequencies obtained by fitting are listed 
in Table 1. We note that, on the one hand, the BiFeOa 
phonons below 150 cm“^ are better resolved in the Eext |1 
[110] polarization with respect to the substrate crys¬ 
tal axes, because the TbScOa phonons are weak in this 
case. On the other hand, the BiFeOa nrodes at 372 and 
446 cm“^ are seen only in the Eext || [001] polarized 
spectra, due to a more favorable TbScOa spectrum in 
this region. Namely, the large TO-LO splitting of the 
TbScOa phonons within 370-480 cin“^ enhances the sen¬ 
sitivity of the IR reflectance to the thin-film phonons— 
Therefore, due to the properties of TbScOa , the two 
polarizations provide complementary information on the 
BiFeOa response. Note that some TbScOa phonons seen 
below 150 cm“^ (see Fig. |3|) exhibit anomalous tempera¬ 
ture shifts and even splitting on cooling. This is probably 
caused by phonon interaction with the crystal field (and 
related electronic transitions) — 

The in-plane lattice parameters of the (110) TbScOa 
substrate are slightly smaller than those of BiFeOa, in¬ 
ducing only a small compressive strain (~ 0.24%) in the 
(OOl)cub epitaxial film. Here, the orientation of the thin 
film is marked with respect to the pseudocubic crystal 
axis (which is denoted by the subscript “cub”) of BiFeOa. 
The ferroelastic and ferroelectric domain structure of 
BiFeOa/TbScOa was investigated in Refs. [Hii where 
two kinds of stripe-like domains separated by (OlO)cub 
vertical boundaries were reported. The spontaneous po- 
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FIG. 3: (Color online) Room-temperature and 10 K IR reflec¬ 
tivity spectra of the TbScOs substrate and reflectance of the 
BiFeOs/TbScOs thin film for polarizations (a,b) Eext || [110] 
and (c,d) Eext jj [001]. 


larization P in adjacent domains is rotated by 109° and 
its direction is tilted from the [001]ciih direction4^ In our 
IR spectra, only vibrations polarized in the film plane 
are active. Since the ferroelectric polarization P is tilted 
from the normal of the thin film plane, phonons of both 
E and Ai symmetries can be excited. Nevertheless, we 
see only some of the modes due to our limited sensitivity 
to the thin film. In contrast to the ceramics samples, the 
phonon damping observed in the thin him is much lower 
(see Table 1). The values of the damping constants are 
comparable to those of single crystals^i which conhrms 
the high quality of our epitaxial thin Him. Similar to the 
case of BiFeOa ceramics, the phonon eigenfrequencies in 
the Him decrease on heating, leading to an increase in the 
static permittivity e(0,T) (see Fig. [2|). The phonon fre¬ 
quencies of the BiFeOa ceramics and the Him present no 
significant differences, which is clearly a consequence of 
the very small strain applied to the film by the substrate. 
The calculated e(0, T) is, however, smaller in the BiFeOa 
film than in the single crystal (see Fig. El), because, for 
the former sample, only strong phonons were revealed in 
the IR reflectance spectra. 

In the low-temperature Eext || [001] spectra, a weak 


but sharp and clearly observable minimum in reflectance 
develops near 600 cm~^ (see Fig. |3][d)). We can ex¬ 
clude its phonon origin, because the factor-group analysis 
in rhombohedral structure does not allow an additional 
mode. A lower crystal symmetry is also excluded be¬ 
cause our thin film has only 0.24% compressive strain; it 
is known that the films change the structure only with 
strain higher than 2%4^ Also the multi-phonon origin 
is not likely since multi-phonon scattering usually de¬ 
cays on cooling; at RT, the peak intensity is markedly 
lower, close to our sensitivity limit (see Fig. |3l[c)). In¬ 
terestingly, the position of the peak, ca. 604 cm“^, cor¬ 
responds to the maximum magnon energy observed out 
of the Brillouin-zone center by inelastic neutron scatter¬ 
ing in a BiFeOs single crystal.— Therefore, we assign 
this peak to scattering by the highest-energy part of the 
magnon branch. The magnon density of states can be ac¬ 
tivated in the IR spectra due to the modulated cycloidal 
magnetic structure of BiFeOa; the maxima in the density 
of states occur either at the maximal energy or at the 
Brillouin-zone boundary. Correspondingly, another peak 
in the density of states can be expected for the magnon 
branch at the K and M points. These magnons have an 
energy^ of ca. 525 cm“^, which falls into the range of 
the highest-energy phonon. We assume that this strong 
phonon screens the weak magnon signal, which is why 
we do not detect the corresponding absorption in the IR 
spectra. Finally, let us mention that the magnons from 
the Brillouin-zone boundary become frequently electri¬ 
cally active due to the exchange striction;^ therefore, 
they can be called electromagnons. 


C. Magnetic excitations in THz spectra of BiFeOs 
ceramics 

Our efforts to measure THz spectra of the BiFeOa film 
failed, due to insufficient signal. This can be explained 
by the fact that the thin film only slightly absorbs the 
THz pulses and the index of refraction n of BiFeOg is 
higher than that of the substrate only by ~ 50%. We 
note that time-domain THz spectroscopy can be used 
only for studies of the thin films with n at least one or¬ 
der of magnitude higher than that of the substrate. As 
an example, such a technique was successfully used for 
investigation of the ferroelectric soft mode behavior near 
a strain-induced ferroelectric phase transition in SrTiOa 
thin fihns4^ 

In contrast, THz spectra of the BiFeOs ceramics were 
successfully measured at temperatures from 5 to 900 K. 
The spectra of the complex refractive index 7V(w) = 
n(uj) -\- ik(uj), determined from experimental data for 
various temperatures, are presented in Fig. |4l In mag¬ 
netic systems, N depends on the complex permittivity e* 
and permeability via the relationship N{u}) = y/fj,*e*. 
Since we cannot resolve whether the observed modes con¬ 
tribute to fi* (as magnons) or e* (as electromagnons or 
polar phonons), we only present n(uj), k(uj) spectra in 
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Wavenumber (cm ’) 

FIG. 4: (Color online) Temperature dependences of the (a) re¬ 
fractive index and (b) extinction coefficient spectra of BiFeOs 
ceramics, determined from the THz transmission. The spec¬ 
tra are shown for every 50 and 100 K below and above RT, 
respectively. 

Fig. m 

In the THz spectra, one can see a gradual increase in n 
and k on heating. This is mainly caused by the softening 
of the TOi phonon and an increase in the phonon damp¬ 
ing with heating. The weak maxima in the k(uj) spectra 
correspond to frequencies of magnetic excitations. Their 
temperature dependence is plotted in Fig.[5]together with 
the modes observed in Raman^SiSi and I R^ 6 , 3 i spectra. 
Note that the modes between 32 and 44 cm“^ become 
IR-active only in external magnetic heldj^i^ Besides a 
good agreement between the frequencies of Raman and 
IR~active modes, one can see a gradual decrease in all 
mode frequencies with increasing temperature. Although 
the excitations should exist at least up to Tn, we fit our 
THz spectra only up to 300 K, because, above RT, their 
damping is very high, precluding their exact fitting. We 
clearly distinguish five magnetic excitations. The three 
of them appearing up to 27 cm”’^ correspond to the IR- 
active modes observed earlier by other authorsji^i 2 £i 2 i 
In our spectra, they exhibit an enhanced damping in 
comparison with single-crystal data, mainly due to the 
fact that we measured un-polarized spectra. At 5 K, 
two modes are seen with peak absorption frequencies of 
53 cm”’^ and 56 cm“^. These modes were not observed in 
IR spectra before, apparently because the samples used in 



Temperature (K) 


FIG. 5: (Golor online) Temperature dependence of the spin- 
excitation frequencies determined by fits of the absorption 
index spectra in the THz range compared with published 
frequencies obtained from Raman scatterin g^”'^^ and far IR 
spectra;^. The modes at ~ 32, 41 and 44 cm“’ were pre¬ 
dicted by Fishmarr^ and they become IR-active only in mag¬ 
netic field above 7T— . 



FIG. 6: (Color online) Magnetic field dependence of the ex¬ 
perimental THz extinction coefficient in the BiFeOs ceramics 
measured at 5K. The dotted line shows the zero-field spec¬ 
trum after applying the magnetic field. 

previous studies were opaque above 40 cm“^. These peak 
frequencies correspond well to the If 7 and $§ modes re¬ 
ported in Raman spectra by Cazayous et al.^ Note that 
Komandin et alr^ predicted a polar and heavily damped 
mode at an estimated frequency of 47cm“^, but their 
rough estimation was based on discrepancies between 
IR reflectivity and THz transmission spectra. Nagel et 
ah^ discovered additional spin excitations near 40 and 
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43 cm“^, but a magnetic field higher than 5 T was needed 
for the activation of these excitations in far-IR transmis¬ 
sion spectra. Furthermore, Fishman predicted^ a mode 
above 50cm“^, which should soften in an external mag¬ 
netic field, and de Sousa et al. predicted^^ exactly the 
two modes which we observe. 

We measured THz transmittance at 5 K under a mag¬ 
netic field up to 7T (see the k{u!) spectra in Fig. [51). 
All modes exhibit small frequency shifts with increas¬ 
ing magnetic field. A similar behavior was reported by 
Nagel et al.^ but their frequency shifts were higher due 
to much higher values of the applied magnetic field (up to 
31 T). We see also weak indications of the modes near 32, 
35, 38, and 42cm~^, which correspond well (except the 
one at 35cm“^) to the magnetic-field induced modes of 
Ref.m In our measurements, while the sample was kept 
at a temperature of 5K, applying the magnetic field of 
7 T irreversibly changed the shape of the magnetic modes 
in the absorption index spectrum (see Fig. [ 5 ]), similarly 
to the observations in Ref. |3l|. This behavior is proba¬ 
bly a consequence of magnetic-field-induced changes in 
the geometry of magnetic domains and pinning the do¬ 
main walls on defects. The pronounced modes near 53 
and 56 cm“^ can be identified with the pair of modes pre¬ 
dicted by Fishman near 45 cm“^ at i7ext = 15Ti^ Thus, 
our observations confirm the complex effective Hamilto¬ 
nian describing the magnetic interactions in BiFeOa; it 
includes nearest and next-nearest-neighbor exchange in¬ 
teractions, two Dzyaloshinskii-Moriya interactions, and 
an easy-axis anisotropyi^ Very recently, this model was 
also confirmed by low-energy inelastic neutron scattering 
spectral^ 

We were not able to fit unambiguously the THz spec¬ 
tra above 300 K using the modes seen at low tempe¬ 
ratures, owing to their heavy damping at high tem¬ 
peratures. Nevertheless, a fit with one effective over¬ 
damped mode provided reasonable results. Surpris¬ 
ingly, this overdamped mode with a relaxation frequency 
around 15cm“^ and Ae = 0.5-1.2 (Ae rising with T) 
was necessary even above Tn, suggesting that there 
are still some paramagnons present in the paramagnetic 
phase. Alternatively, this mode might be due to a multi¬ 
phonon or quasi-Debye absorption allowed in the non- 
centrosymmetric phase4^ 

Finally, we would like to stress that from unpolarized 
THz spectra, we could not distinguish whether our mag¬ 
netic excitations are magnons or electromagnons. The 
symmetry analysis published by Szaller et al.^ how¬ 
ever, allows electromagnons in the cycloidal G-type an¬ 
tiferromagnetic phase of BiFeOa. In analogy with selec¬ 
tion rules for polar phonons which, in acentric lattices, 
are both IR and Raman active, we suggest that for spin 


excitations, their simultaneous IR and Raman activites 
indicate that they are electromagnons. During the re¬ 
view process, the polar activity of spin excitations has 
been confirmed by Kezsmarki et al.^ who observed di¬ 
rectional dichrois m^^i^^ of a BiFeOa single crystal in po¬ 
larized far-IR spectra between 10 and 30cm“^. 

IV. CONCLUSIONS 

An extensive study of IR vibrational spectra of BiFeOa 
ceramics and an epitaxial thin film is reported. The 
intensities of all phonons observed in the ceramics are 
higher than in the previous publications Thus, the 
static permittivity of our samples calculated from phonon 
contributions is close to the previously published single 
crystal dataJi Some phonons slightly soften on heat¬ 
ing, leading to an increase in the permittivity towards 
Tq. Nevertheless, the permittivity is much lower than in 
canonical ferroelectric perovskites, such as BaTiOa and 
KNbOa, because the phonons in BiFeOa are much harder 
and the Born effective charges are much smaller. In ad¬ 
dition, the phonons in a BiFeOa/TbScOa epitaxial thin 
film were studied for the first time, showing parameters 
similar to those in BiFeOa single crystals. In the thin 
film, an additional weak excitation near 600 cm~^ was 
detected which apparently corresponds to a peak in the 
magnon density of states and which can be excited due to 
the incommensurately modulated magnetic structure of 
BiFeOa. Time-domain THz spectra of BiFeOa ceramics 
reveal most of the spin excitations previously observed in 
single crystalsj ^^i^^ Also, at 5 K, a pair of IR-active modes 
near 55cm“^ were observed, which corresponds to spin 
excitations theoretically predicted by de Sousa et alM- 
and Fishmani^ This observation confirms the particular 
form of the Hamiltonian suggested by Fishman for the 
explanation of the magnetic interactions in BiFeOa. 
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